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Abstract 



We demonstrate the coupling of bright and dark Surface Lattice Resonances (SLRs), which are 
. collective Fano resonances in 2D plasmonic crystals. As a result of this coupling, a frequency stop- 



gap in the dispersion relation of SLRs is observed. The different field symmetries of the low and 
high frequency SLR bands lead to pronounced differences in their coupling to free space radiation. 
Standing waves of very narrow spectral width compared to localized surface plasmon resonances 



are formed at the high frequency band edge, while subradiant damping onsets at the low frequency 
band edge leading the resonance into darkness. We introduce a coupled oscillator analog to the 
plasmonic crystal, which serves to elucidate the physics of the coupled plasmonic resonances and 
to estimate very high quality factors (Q > 700) for SLRs, which are the highest known for any 2D 
plasmonic crystal. 
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Metallic nanoparticles supporting surface plasmon resonances allow light to be localized 
in nanoscale volumes, thereby opening exciting possibilities such as nanoscale control of 
emitters large electromagnetic enhancements {2], and nonlinear nano-optics J3]. Much 
attention has been given to Localized Surface Plasmon Resonances (LSPRs), which arise 
in individual particles when their conduction electrons are coherently driven by an electro- 
magnetic field. Although localized surface plasmons may couple, their resonances are in 
general severely broadened due to strong radiative damping and hence exhibit low quality 
factors Q. A recent development in nanoplasmonics deals with collective resonances in pe- 
riodic arrays of metallic nanostructures, or plasmonic crystals. Such arrays support Surface 
Lattice Resonances (SLRs), which are collective resonances mediated by diffractive coupling 
of localized plasmons. This coupling occurs near the critical frequency when a radiating 
diffraction order becomes evanescent, i.e., at the Rayleigh anomaly. SLRs were introduced 
by Carron (4], and the interest in this phenomenon was revived by Schatz and co-workers 
with a series of works on ID and 2D arrays 5|, |6[]. However, the experimental observation 
of SLRs was elusive for many years {tJ • Recent advances in nano-fabrication and in the un- 



lowed for their observation in periodic arrays of nanostructures 



derstanding of SLRs have a 

with different geometries fe. In contrast with LSPRs, SLRs possess mueh h,ghe r Qs »d 
the associated polaritons can propagate over tens of unit cells in the plasmonic crystal |12j . 
The relevance of SLRs for enhanced, directional, and polarized light emission and 
sensing 15[ has been recently demonstrated. Although the coupling of surface modes in 
periodic metallic structures has attracted much interest 16l-ll9|. especially for its connection 
with frequency stop-gaps [20], coupled SLRs have not been discussed yet. 

In this paper, we demonstrate the mutual coupling of SLRs and the formation of a 
frequency stop-gap in the dispersion relation of these modes. This coupling leads to a strong 
modification of the SLRs characteristics, including the onset of subradiant damping in the 
low frequency band, zero group velocity modes in the high frequency band, and Q-factors 
for both bands which are amongst the highest reported for any 2D plasmonic crystal. Our 
results set the basis for controlling the dispersion of SLRs, and they open new possibilities 
in sensing, enhanced spontaneous light emission, and lasing at the band edges of SLR gaps. 

We have investigated 3x3 mm 2 arrays of gold nanorods fabricated on a silica substrate 
using substrate conformal imprint lithography (SCIL) 2l|. A top view SEM image of an 
array is displayed in the inset of Figured] This array has rods with dimensions 450 x 120 x 38 



nm 3 arranged in a lattice with constants a x = 600 nm and a y = 300 nm. The array was 
embedded in a uniform surrounding medium by placing a silica superstrate preceded with 
n=1.45 index matching fluid to ensure good optical contact. We measured the variable angle 
transmittance of the collimated beam from a halogen lamp while rotating the sample around 
the y-axis. The polarization of the incident light was set along the y-axis, probing the short 
axis of the nanorods. For this polarization the dipolar LSPR lies at higher energies than 
the (±1,0) diffraction orders, thereby allowing the coupling of localized surface plasmons to 
these orders 

Figure [T] displays the extinction of the array defined as 1 — T, with T the transmittance, 
as a function of the reduced frequency, i.e., the angular frequency normalized by the speed 
of light in vacuum, and the projection of the incident wave vector onto the surface of the 
array k\\ = - sin(#)x, with 9 the angle of incidence. The broad, dispersionless extinction 
peak centered at a frequency near 9 mrad/nm corresponds to the excitation of LSPRs in 
the individual nanorods. The sharp, dispersive peaks at lower frequencies correspond to the 
excitation of SLRs. The (±1, 0) Rayleigh anomalies, which are the conditions for which the 
(±1, 0) diffracted orders are grazing to the surface, are indicated with solid lines in Figured] 
The coupling of localized surface plasmons to the Rayleigh anomalies is the origin of the 
observed SLRs. 

A salient feature in the measurements of Figured] is the formation of a stop-gap centered 
at u/c = 6.85 mrad/nm and near k\\ = 0, where the two SLRs mutually couple. We note 
that this is not a complete photonic bandgap, since it exists for y-polarized light only. The 
gap arises from the coupling of two counter-propagating surface polaritons which, due to 
the structural anisotropy of both the nanorods and the lattice, have a strong polarization 
dependence on their coupling to free space radiation. At the high frequency band edge, we 
observe that the dispersion of the (+1,0) SLR flattens. This flattening of the band can be 
translated as a reduction of the mode's group velocity and the formation of standing waves, 
which are also associated with an increased density of optical states. At the low frequency 
band edge, the (-1,0) SLR becomes weaker and narrower. This behavior is characteristic of 
a mode tending towards subradiance, where radiative damping is suppressed in a collective 



state with an antisymmetric wave function 22]. As shown in nanoslit arrays 18[, there is an 



intimate connection between subradiant damping and the opening of a gap in the dispersion 
relation of surface plasmon polaritons (SPPs). 



Figure [2](a) shows a close view of the stop-gap in Figure [TJ and Figure 2(b) shows re- 
sults from finite element simulations (COMSOL). For the simulations we used a constant 



refractive index of 1.45 for silica, the permittivity of gold as given in Ref. [23|, and Bloch- 
Floquet boundary conditions. The transmittance was calculated as the ratio between the 
transmitted power through the array and the incident power. In Figures [2£c) and (d) we 
compare simulations with measurements at k\\ = and k\\ = 0.4 mrad/nm, respectively. By 
reproducing the measured dispersion of SLRs and the gap's central frequency and width, a 
good qualitative agreement between measurements and simulations is demonstrated. Dis- 
crepancies in the amplitude and spectral width of the resonances can be mainly attributed 
to differences between the simulated and fabricated geometries, especially near the corners. 

It can be appreciated in Figures [1] and |2] that the SLR peaks together with the Rayleigh 
anomaly dips give rise to asymmetric resonance lineshapes, which can be understood in the 



framework set forth by Fano |24j. Fano described the quantum interference between a dis- 
crete state and a continuum of states as the origin of asymmetric resonance lineshapes. It 
was later shown in metallic subwavelength hole arrays that the coupling of surface plasmons 
to Rayleigh anomalies leads to similar lineshapes, with the resonantly scattered light acting 



as the discrete state and the background transmission as the continuum [25J. A similar sit- 
uation is observed in our configuration. The broad dipolar LSPR determines the extent of 
background transmission, i.e. the continuum, according to its frequency difference with the 
Rayleigh anomaly. The Rayleigh anomaly corresponds to the discrete state. The diffractive 
coupling of localized surface plasmons therefore resembles the interaction between a contin- 
uum of states and a discrete state, leading to asymmetric resonance line shapes. As seen 
in Figures [1] and [21 the degree of asymmetry of these line shapes changes depending on the 
frequency and k-vector of excitation. This dependency is rooted in the relative contributions 



of resonant and non resonant scattering |26|], which also manifest as a modification of the 
SLR damping; the latter point will be addressed further in the text. 

The different electrodynamic response leading to the bright and dark character - efficient 
and inefficient coupling to light - of the (+1,0) and (-1,0) SLRs, respectively, transpires from 
the near-field enhancement and surface charge distribution of the nanorods at the respective 
frequencies. In Figure |3] we show simulation results for ku = 0.4 mrad/nm at two frequencies: 
u/c = 7.1 mrad/nm and oj/c = 6.7 mrad/nm, which correspond to the (+1,0) and (-1,0) 
SLR, respectively. The small angle of incidence was chosen such that the extinction is not 



negligible for the (-1,0) SLR. Both plots are at a plane parallel to the array located at the 
mid-height of the nanorods. Charges of opposite sign at the surface of the nanorods are 
plotted in white and black, while the total near-field enhancement, i.e., | _£7 1 2 / 1 -E'o 1 2 with E 
the total field and Eq the incident field, is displayed by the color scale. Figure [3] illustrates 
how the different resonant response has its origin in the symmetry of the modes. In order 
to couple to the incident plane wave at normal incidence, the mode has to be symmetric 
with respect to the plane defined by the incident k- and polarization vectors intersecting 
the nanorods along their center, i.e. the symmetry plane indicated by the dotted lines in 
Figure El The (+1,0) mode, which is shown in Figure El^a) for a small angle of incidence, 
has symmetric field and charge distributions with respect to the symmetry plane. A strong 
dipole moment is seen for each nanorod, and a strong dipolar inter-rod coupling takes places 
along the y-direction also. This results in a large extinction and an efficient coupling of the 
(+1,0) mode to normal incidence light, as it can be appreciated in Figure [TJ In contrast, the 
(-1,0) mode has an antisymmetric field and charge distribution for k» = 0; the net dipole 
moment is therefore zero and the extinction vanishes. This symmetry is broken for angles of 
incidence larger than 9 = 0°, thus allowing the excitation of this resonant mode as shown in 
Figure [3](b) . The broken symmetry manifests as a quadrupolar surface charge distribution 
displaced from the symmetry axis. This results in a nonzero intra-rod and inter-rod dipole 
moment, which can be recognized from the charges of opposite sign inside the nanorods and 
for adjacent nanorods along the symmetry axis, respectively. 

The coupled nature of SLRs can be elucidated by making an analogy with a set of three 
mutually coupled harmonic oscillators. Coupled oscillators have proven useful in under- 



standing electromagnetic phenomena 27], |28[. In this analogy, the conduction electrons in 



the nanorods driven by the electromagnetic field are modeled as oscillator 1 driven by a 
harmonic force F = F e~ lulBt , whereas the (+1,0) and (-1,0) Rayleigh anomalies are modeled 
by oscillators 2 and 3, respectively. The equations of motion for the system are 
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where Xj, jj, and Uj (j = 1,2,3) are the the displacement from equilibrium posi- 
tion, damping, and eigenfrequency associated with the j th oscillator, respectively, and fL fc 



TABLE I. The Qjk terms (j, k = 1, 2, 3 and j ^ k) are the coupling frequencies between the j and 
k th oscillators, and the 7j terms are the damping frequencies associated with the j th oscillator, for 
the system described by Equation 1. All quantities are given in units of mrad/nm. In the entries 
for which a minimum estimate is given, the value in parenthesis represents the value yielding the 
spectra in Figure 4. 
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=0 


2.8 ± 0.1 






<0.01 (0.001) 




k\ 


=0.17 


3.1 ± 0.1 


<0.6 (0.1) 


1.1 ± 0.1 


0.02 ± 0.01 


0.020 ± 0.005 


k\ 


=0.68 


3.3 ± 0.1 


2.4 ± 0.1 


<0.7 (0.3) 


0.06 ± 0.02 


0.03 ± 0.01 



{k = 1,2,3 and j ^ k) is the coupling frequency between the j th and k th oscillator. Since 
we are interested in the extinct optical power in driving the electrons in the nanorod, we 
calculate the absorbed mechanical power by oscillator 1 from the driving force, which is 
given by P(t) = Fx\. Integrating P(t) over one period of oscillation and scanning the driv- 
ing frequency u s yields an absorbed power spectrum P(oj s ), which is representative of the 
extinction spectrum. 

In Figure 4 we compare P(u s ) with the measured extinction spectra at three values of 
k\\. In all three cases oj\ = 9.25 mrad/nm and 71 = 2.3 mrad/nm, which reproduce the 
frequency and damping of the LSPR, and Fq = 0.57 units of force per mass is a fitting 
parameter determining the amplitude of the spectra. The eigenfrequencies U2 and U3 are 
determined for each value of k\\ by the corresponding (±1,0) Rayleigh anomalies of the 
array. The remaining parameters, i.e., coupling and damping frequencies associated with 
the lattice modes, are the parameters used to fit the measured SLR lineshapes and are given 
in Table 1 for three values of k\\ . 

Figure H^a) displays the spectra at k» = 0, where the Rayleigh anomalies are degenerate 
at U2 = 7.26 mrad/nm and the (-1,0) SLR is a dark state. In this case the model reduces 
to that of two coupled oscillators. With increasing ku the (-1,0) SLR comes out of the 
darkness, so the three oscillators are mutually coupled. In Figure S(b) we consider the 
case k\\ = 0.17 mrad/nm, where U2 = 7.33 mrad/nm and U3 = 6.92 mrad/nm. From the 
values given in Table 1 we see that with respect to the normal incidence case, at k» = 0.17 
mrad/nm the nanorods are more strongly coupled to the (+1,0) Rayleigh anomaly (f2i2 
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increases), only weakly coupled to the (-1,0) Rayleigh anomaly (low ^13), the Rayleigh 
anomalies are mutually coupled (high Q23), and the damping of both resonances has very 
significantly increased (both 72 and 73 increase). Further increasing to kn = 0.68 mrad/nm 
makes U2 = 7.68 mrad/nm and W3 = 6.70 mrad/nm, which is the case in Figure 4(c). At 
this value of kn there is an increased coupling of the nanorods to both Rayleigh anomalies 
(both Q12 and O13 increase), the coupling between the Rayleigh anomalies decreases 
decreases), and the damping of both resonances increases (both 72 and 73 increase). 

From the above mentioned behaviors and the quantities given in Table 1 we are able to 
draw several conclusions. Firstly, the model shows how the coupling terms Vtjj, determine 
the frequency difference between the SLR peaks and the Rayleigh anomalies as evidenced in 
Figure [TJ Namely, for the (+1,0) Rayleigh anomaly we observe that as k\\ increases the SLR 
peak deviates more in frequency. This behavior is defined in the increase of Q12 which detunes 
the resonance peak from the Rayleigh anomaly. In contrast, for the (-1,0) diffraction order 
we observe a decreasing frequency deviation of the SLR from its corresponding Rayleigh 
anomaly as kn increases. In this case the dominant interaction near normal incidence is the 
mutual coupling of SLRs described by the term O23, which detunes the (-1,0) SLR from 
its Rayleigh anomaly at low k\\. Although Q 13 > Q 2 3 at large values of kn, at low values 
of k\\ the Q23 interaction dominates due to the smaller frequency difference between the 
eigenfrequencies U3 and 002- Secondly, the model shows that the damping of both resonances 
increases with kn , which leads to less asymmetric Fano lineshapes and broader linewidths for 
both resonances. The decreasingly asymmetric lineshapes as kn increases are especially clear 
in Figure HJ whereas the SLR broadening and variable extinction of the Rayleigh anomalies 
are visible in Figure [TJ 

It is interesting to calculate the Q-factors of the uncoupled oscillators, which follow 
from the definition Qj = 0Uj/jj. For oscillator 1 resonating at the LSPR frequency we 
obtain Qi = 4. For oscillator 2 yielding the (+1,0) SLR, we have Q2 > 700 at normal 
incidence. This value, which is to the best of our knowledge higher than any reported Q 
for an experimental study on 2D plasmonic crystals, is a minimum limit. This limit arises 
because for any value 72 < 0.01 a fit to the measurement within a 10% uncertainty in the 
magnitude of the extinction can be obtained. For oscillator 3 yielding the (-1,0) SLR, we 
have Q3 = 300 at kn = 0.17 mrad/nm. It is important to realize that the resonances in 
the coupled system exhibit an effective damping which is not equal to the damping of the 
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uncoupled oscillators. Nevertheless, the above values clearly reflect the large differences in 
Q-factors between LSPRs and SLRs. Furthermore, the coupled oscillator model points to the 
origin of the narrow SLR linewidths, which is the coupling of two harmonic oscillators with 
very different damping. By comparing Qi and Q3 with previously reported experimental 
data, we get an insight into how well the Q-factors of the uncoupled oscillators represent 
the Q-factors of the resonances in the coupled system. Firstly, Qi = 4 is a typical result 



for LSPRs with high radiative damping 



Ropers and co-workers 



29] . Secondly, in Ref. 

measured for subradiantly damped SPPs in a nanoslit array near a stop-gap (similar to 
what we observe for the (-1,0) SLR near the gap) lifetimes near r = 200 fs. Nanoslits can 
be represented quite accurately as resonators |30j, so the lifetime of the excited state, r, is 
related to Q by Q = 2nur, with v the frequency Inserting r = 200 fs and v = 400 THz (the 



frequency of the subradiant mode in Ref. 18J) yields Q = 500, which is comparable with the 
value that we find, i.e., Q3 = 300. Despite the differences in the structures considered and 
possibly the extent of subradiant damping, it is remarkable that the oscillator model yields 
at the very least an order of magnitude estimate for the Q of a nearly subradiant SPP mode 
in the vicinity of a stop-gap. It should be mentioned that for arrays of dimensions on the 
order of 10 microns or less, the Q-factors of SLRs are expected to appreciably decrease. This 
length scale corresponds to the associated surface polariton propagation lengths, which were 



recently determined in Ref. [12|. For arrays much larger than this characteristic length, such 
as the one herein considered, the Q-factor saturates. The inhibition of radiative damping 
that leads to a narrowing of the SLR linewidth is therefore a collective effect, i.e. it depends 
on the number of particles present in the array. 

In conclusion, we have demonstrated the coupling of SLRs and the associated opening 
of a frequency stop-gap in the dispersion relation of these modes. The symmetric and 
antisymmetric field/charge distributions responsible for the bright and dark nature of the 
(+1,0) and (-1,0) modes, respectively, were illustrated. We have also estimated the very 
high quality factors of SLRs (Q > 700), which are to the best of our knowledge the highest 
reported values for any experimental study on 2D plasmonic crystals. Coupled bright and 
dark collective modes, as well as stop-gaps, offer the possibility to carefully design abrupt 
changes in the Local Density of Optical States (LDOS) over narrow spectral regions and in 
an extended volume. LDOS manipulation in plasmonic structures is relevant for enhanc- 
ing the efficiency of light emitting devices (LEDs), sensors, and nonlinear processes, all of 
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which can be tailored in a frequency, angle, and polarization dependent manner with cou- 
pled SLRs. Moreover, the standing waves formed at the high frequency band edge hold 
exciting properties for the manipulation of light at the nanoscale since they have zero group 
velocity while being much less damped than LSPRs. Finally, we envisage that the strong 
suppression of radiative losses herein discussed holds great promise for the development of 
high Q distributed feedback surface polariton lasers and plasmonic sensors with enhanced 
sensitivity. 
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k\\ (mrad/nm) 

FIG. 1. (Color) Extinction spectra as a function of kn for y-polarized light incident on the array 
of gold nanorods shown in the inset. The black lines with positive and negative slope indicate the 
(+1,0) and (—1,0) Rayleigh anomalies, respectively. The dispersionless resonance at 9 mrad/nm 
is the dipolar localized surface plasmon resonance for the short- axis of the nanorods, whereas the 
narrower and dispersive resonances below the Rayleigh anomalies are the surface lattice resonances. 
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FIG. 2. (Color) Measurements (a) and finite element simulations (b) of the extinction spectra 
of the array shown in Figure 1. Figures (c) and (d) are cuts at kn = and kn = 0.4 mrad/nm, 
respectively, of both (a) and (b) . The open circles in (c) and (d) are measurements and the solid 
curves are simulations. 
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FIG. 3. (Color) Near field enhancement in color scale and surface charge distribution (at an 
arbitrary phase) in black and white at the mid-height of the nanorods for the (+1,0) (a) and (-1,0) 
(b) surface lattice resonance at h\ = 0.4 mrad/nm. The dotted lines indicate the plane of symmetry 
for coupling to radiation. 
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FIG. 4. The black open circles are cuts of the measured extinction spectra shown in Figure 1 at 
three values of fen: (a) fcy = mrad/nm, (b) k\\ = 0.17 mrad/nm , and (c) k\\ = 0.68 mrad/nm. 
The blue solid curves represent the absorbed power in oscillator 1 of the coupled oscillator model 
described in the text, with coupling and damping frequencies as given in Table 1. 
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